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Abstract 

The Dalitz plot of rj' r]Tr~^7T~ decay is studied using (225.2 ± 2.8) X lO'^ J/ip events collected 
with the BESIII detector at the BEPCII e'^e~ collider. With the largest sample of rj' decays to 
date, the parameters of the Dalitz plot are determined in a generalized and a linear representation. 
Also the branching fraction of J/i/j 777' is determined to be (4.84 it 0.03 it 0.24) x 10~^, where 
the first error is statistical and the second systematic. 

PACS numbers: 12.39.-x, 13.25.Gv, 14.40.Be 
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I. INTRODUCTION 



Chiral Perturbation Theory is the low energy effective theory of Quantum Chromody- 
namics. Below the p mass region, the interactions of the (tt, K, 77) particles are systematically 
analyzed within this framework. The success in the description of these low-energy inter- 
actions makes ChPT a powerful theoretical tool [l]. Although the mass of the r/' is high 
and 7]' — i- 'r]TT~^TT~ decay has a low Q value, which limit the predictive power of the Effective 
Chiral Lagrangian model, the experimental study of the process may supply information to 
test the predictions of chiral theory 2|- 4 ] and possible extensions of ChPT such as large-NC 
ChPT and resonance Chiral Theory 5j. The hadronic decays of the rj' meson have also 
been extremely valuable in studies devoted to the effect of the gluon component [G] and 
the possible nonet of light scalars [7]. Previously, the GAMS-47r and VES Collaborations 
have measured the related Dalitz plot parameters (GAMS-47r for the r]' — )■ r^vr^TT^ channel [8| 
and VES for rj' — > r^vr+vr" {q]) complementing older results reported by an early GAMS [lO| 
and CLEO [11] Collaborations. In the isospin limit, the values of the Dalitz plot parame- 
ters should be the same; however the experimental measurements show some discrepancies 
among them. 

In this article, with a new level of precision, we present results for the Dalitz plot pa- 
rameters for 1]' —7- i]7i^7i^ based on (225.2 ± 2.8) x 10^ J/ip events collected by BESIII at 
BEPCII. 



II. BESIII AND BEPCII 



BESIII/BEPCII [12^ is a major upgrade of the BESII experiment at the BEPC accel- 
erator |l^ for studies of hadron spectroscopy and r-charm physics The design peak 
luminosity of the double-ring e~^e~ collider, BEPCII, is 10^^ cm"^s~^ at a beam current of 
0.93 A. The BESIII detector with a geometrical acceptance of 93% of Att, consists of the fol- 
lowing main components: 1) a small-celled, helium-based main draft chamber (MDC) with 
43 layers. The average single wire resolution is 135 m, and the momentum resolution for 1 
GeV/c charged particles in a 1 T magnetic field is 0.5%; 2) an electromagnetic calorimeter 
(EMC) made of G240 Csl (Tl) crystals arranged in a cylindrical shape (barrel) plus two 
endcaps. For 1.0 GeV photons, the energy resolution is 2.5% in the barrel and 5% in the 



endcaps, and the position resolution is 6 mm in the barrel and 9 mm in the endcaps; 3) a 
Time-Of-Flight system (TOF) for particle identification composed of a barrel part made of 
two layers with 88 pieces of 5 cm thick, 2.4 m long plastic scintillators in each layer, and 
two endcaps with 96 fan-shaped, 5 cm thick, plastic scintillators in each endcap. The time 
resolution is 80 ps in the barrel, and 110 ps in the endcaps, corresponding to better than 
a 2 sigma K/n separation for momenta below about 1 GeV/c; 4) a muon chamber system 
(MUG) made of 1000 m^ of Resistive Plate Chambers (RFC) arranged in 9 layers in the 
barrel and 8 layers in the endcaps and incorporated in the return iron of the superconducting 
magnet. The position resolution is about 2 cm. 

The estimation of physics backgrounds are performed through Monte Carlo (MC) simu- 
lations. The GEANT4-based simulation software BOOST |l^ includes the geometric and 
material description of the BESIII detectors, detector response and digitization models, as 
well as the tracking of the detector running conditions and performance. The production 
of the J /ip resonance is simulated by the MC event generator KKMC IGj, while the decays 
are generated by EvtGen [l7| for known decay modes with branching fractions being set 
to the FDG [18] world average values, and by Lundcharm 1^ for the remaining unknown 
decays. The analysis is performed in the framework of the BESIII Offline Software System 
(BOSS) [20] which takes care of the detector calibration, event reconstruction and data 
storage. 



III. EVENT SELECTION 

The 7]' is identified by its decay into tjit^tt^ with 77 —t- 77 in J/-?/; radiative decays, 
and candidate events with the topology 7777r+7r^ are selected using the following criteria. 
Charged tracks in BESIII are reconstructed from MDC hits. To optimize the momentum 
measurement, we select tracks in the polar angle range |cos^| < 0.93 and require that 
they pass within ±10 cm of the interaction point in the beam direction and within ±1 
cm in the plane perpendicular to the beam. Electromagnetic showers are reconstructed by 
clustering EMC crystal energies. Efficiency and energy resolution are improved by including 
energy deposits in nearby TOF counters. Showers identified as photon candidates must 
satisfy fiducial and shower-quality requirements. The minimum energy is 25 MeV for barrel 
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showers (|cos0| < 0.8) and 50 MeV for endcap showers (0.86 < |cos6'| < 0.92). Photons 
in the region between the barrel and endcaps are not well measured and are not used. To 
exclude showers from charged particles, a photon must be separated by at least 20° from 
any charged track. EMC cluster timing requirements suppress electronic noise and energy 
deposits unrelated to the event. 

The TOF (both Endcap and Barrel) and dE/dx measurements for each charged track 
are used to calculate Xp/z)(^) values and the corresponding confidence levels Prohpioii) for 
the hypotheses that a track is a pion, kaon, or proton, where i {i = 'k/K/p) is the particle 
type. For pion candidates, we require ProhpiDij^) > ProhpiD^K) and Prohpioijx) > 0.001. 

Candidate events must have two charged tracks with zero net charge, and the number of 
photons should be greater than two. At least one charged track must be identified as a pion. 
We do four-constraint (4C) kinematic fits imposing energy and momentum conservation 
under the J/ip ^ 7777r+7r^ hypothesis looping over all photon candidates, and select the 
combination with the minimum x^(7777r+7r^). The minimum x^(7777r"'"7r~) should be less 
than 200, and the efficiency of this requirement is around 99%. The r] candidates are selected 
from the combination with the two photons' invariant mass closest to t] nominal mass. With 
the above event selections, a very clear r] signal is observed. In the analysis below, we define 
the ?7 signal region as 0.518 < m^^ < 0.578 GeV/c^, and the rj mass sidebands region as 
0.443 < m^^ < 0.473 GeV/c^ or 0.623 < m^^ < 0.653 GeV/c^. 

The backgrounds in the selected event sample from a number of potential background 



channels listed in the PDG [18| are studied with MC simulations. The background level is 
very low in the rj' mass region. The main backgrounds are from J/tp — > 77p'^ — >■ 
777r+7r~ and Jj^l) — )■ 777' — )■ 770; — > 777r+7r^7r°, which can be described by the nor- 
malized ri mass sidebands events. The other backgrounds with ri candidates are from 
Jj^ 7/i(1285)/?7(1405/1475)//i(1510) -> 7r77r+7r". None of these backgrounds give 
peaking backgrounds in the 77' mass region. The total background contamination is esti- 
mated to be only 0.57% within the Vj mass region (~3cr). An inclusive MC event sample is 
also used to investigate other possible surviving background events, but no other possible 
background from the inclusive MC is found. 



7 



IV. NUMBER OF J/ip EVENTS 



The number of J /if) events, Nj/^, used in this analysis is determined from the number 
of inclusive events. Charged tracks are selected requiring their points of closest approach to 
the beam line be within 15 cm of the interaction point along the beam line and within 1 cm 
in the plane perpendicular to the beam line, their angles with respect to the beam line, 6, 
must satisfy | cos6\ < 0.93, and their momenta must be less than 2.0 GeV/c. Clusters in 
the EMC must have at least 25 (50) MeV of energy in the barrel (endcap) EMC, and have 
Icos^l < 0.93. 

Event selection requires at least two charged tracks and visible energy, Eyis, greater than 
1.0 GeV. Here E^s is defined as the sum of charged particle energies computed from the track 
momenta assuming pion masses, plus the neutral shower energies measured in the EMC. To 
reduce backgrounds from Bhabha and dimuons, events with only two charged tracks must 
have the momenta of the charged tracks less than 1.5 GeV/c and their energy deposit 
in the EMC less than 1.0 GeV. Backgrounds from Bhabha and dimuon events surviving 
the selection criteria are small. The continuum contribution (e+e^ — t- anything) and the 
surviving backgrounds are removed by subtracting the number of events selected with the 
above criteria from a continuum sample taken at a center of mass energy of 3.08 GeV and 
normalized by relative luminosity and the cross section assuming a 1/s dependence. 

The number of J/ if) inclusive events is also determined from the distribution of which 
is the average of the z distances from the interaction point along the beam of the point of 
closest approach of tracks to the beam line. Here the number of J /if) inclusive events is 
taken to be the number of events in a signal region (—4 < z < A cm) minus the number of 
events in sideband regions (6 < |^| < 10 cm) of z. 

The efficiency is determined from data using J /ip events from if)' — )■ TT^Ti~ J /ip decays 21| 
in the BESIII 106 M ip' sample 22|]. MC simulation is used to determine a small correction 
(1.0108) to this efficiency arising from the two extra tracks and the motion of the J/ip 
in the ip' events. This procedure is less sensitive to differences between data and MC 
simulation than using only MC to determine the efficiency. The agreement between data 
and MC simulation is shown for the cos 9 distribution of charged tracks in Fig. [U^a) and the 
total energy deposit in the EMC, Eemc, in Fig- Wljo)- The discrepancy between data and 
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MC simulations in Fig. [U^b) is due to the imperfect MC generator and imperfect detector 
simulation. The systematic error due to the i?,,,., requirement is negligible. 




cose Emc (GeV) 

FIG. 1: (a) The cos^ distribution of charged tracks for events satisfying selection criteria, (b) The 
distributions of the total energy in the EMC for events satisfying selection criteria. Dots are data, 
and the histogram is J/-^ — t- inclusive simulated events. 

The result is Nj/^ = (225.2 ± 2.8) x 10^, where the error is systematic and is determined 
mostly by the track efficiency difference between data and MC (0.41%), the variation with 
the minimum charged track multiplicity requirement (0.78%), the difference when the noise 
levels in the two samples of J/tp and ip' events are modified (0.49%), the error associated 
with fitting the distribution of mass recoiling from the tt+tt^ to determine the number of 
ip' — )■ 71^71^ J /ip events (0.45%), the error due to the continuum subtraction (0.18%), the 
difference between the continuum subtraction and the sideband subtraction methods for 
determining the number of events (0.18%), and the difference for changing the generator 
(0.49%). The statistical error is negligible. A second analysis determines Nj/^ from J/ip 
events, where / is a or e, and obtains consistent results. 



V. BRANCHING FRACTION MEASUREMENT 

Figure [2] shows the invariant mass distribution of r^Tr+vr" candidate events. This distribu- 
tion is fitted with a double- Gaussian function for the rj' signal and a linear function for the 
background shape. The fit yields 43826 ±211 events. The J /ip 777' branching fraction is 
calculated using 



B{J/^ ^ 7^') 



Nji^ X e X B{i]' — )■ r^TT+TT^) x B{i] — )■ 77) ' 



where A^°^* is the number of events observed, Njj^ is the number of J /if) events, and e is the 
selection efficiency obtained from MC simulation, which is 23.57%. The branching fraction is 
then determined to be (4.84 ±0.03) x 10^^, where the error is statistical only. We also check 
B{J/il) — 7- 77]') by using the number of events after the 6C kinematic fit requirement (the 
reconstructed momenta of two gammas are constrained to the r] mass and the reconstructed 
momenta of rjn^Ti^ is constrained to the rj' mass), where the number of signal events is 
obtained by subtracting all the simulated normalized backgrounds with rj candidates and 
normalized rj mass sidebands events directly. The difference for B{J /ip — )■ 777') is only 0.3%. 

10000 
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M(ti tC n ) (GeV/c^) 

FIG. 2: The r]'ir~^7r~ invariant mass distribution of the final candidate events. The dots with error 
bars represent data, and the solid curve is the result of the fit described in the text. The dashed 
curve is the background polynomial. 




VI. MEASUREMENT OF THE MATRIX ELEMENT 

The internal dynamics of the decay rj' — )■ rj-K^-K' can be described by two degrees of 
freedom since all the particles are spin zero particles. The Dalitz plot distribution for the 
charged decay channel rj' — )■ rj-n-^n^ is described by the following two variables: 

where Tjr^r; denote the kinetic energies of mesons in the t]' rest frame and Q = Ty^+TT^+ +Tt^- = 
rrin' — rrij^ — 2171^,. The squared absolute value of the decay amplitude is expanded around 
the center of the corresponding Dalitz plot in order to obtain the Dalitz slope parameters: 

= A{l + aY + bY^ + cX + dX^), (2) 
10 



where a, b, c and d are real parameters and A is a normalization factor. This parametrization 
is called the general decomposition. The parametrization in Eq. ([2]) has also been proposed 
with an extra term, either eXY or fX^ + gY^. For the charged channel t]' t^tt+tt^, odd 
terms in X are forbidden due to charge conjugation symmetry, while for the neutral channel 
7]' — ^ ?77r°7r°, c = from symmetry of the wave function. The Dalitz plot parameters may 
not be necessarily the same for charged and neutral decay channels. However, in the isospin 
limit they should be the same. 

A second parametrization is the linear one [18|: 



= A{\l + aY\^ + cX + dX^), (3) 

where a is a complex parameter. Of particular interest is the real component of the complex 
constant a, which is a linear function of the kinetic energy of the r]. A non-zero value of a 
may represent the contribution of a gluon component in the wave function of the rj' in the 
dynamics of its decay 1^. Comparison with the general parametrization gives a = 2Re(a) 
and b = Re^(a) + Im^(a). Both parameterizations are equivalent if 6 > 

To improve the rj and r]' mass resolutions and reduce the migration of events to the nearby 
bins in the Dalitz plot, we use kinematic information after a 6C kinematic fit to calculate 
the X and Y values. 

Figure [3] (a) shows the experimental form of the Dalitz diagram for the decay r]' — riTT^TT~ 
in terms of the variables X and Y with the rj-K^-K^ mass in the 0.93-0.98 GeV/c^ mass region, 
while the corresponding projections on variables X and Y are shown in Figs. [3] (b) and (c), 
respectively. In Figs. [3] (b) and (c), the dashed histograms are from MC signal sample 
with 7]' — > rj-K^-K^ events produced with phase space, while the solid histograms are the 
fitted results described below. The resolutions in the variables X and Y over the entire 6C 
kinematical region are ax = 0.03 and ay = 0.025, respectively, according to MC simulation. 

The dependence of the matrix element on each variable, X and Y , after integration over 
the other, and after dividing by phase space, is shown in Fig.|H Fitting the data with Eq. ([3]) 
gives the following values of the parameters: Re(a) = —0.035 ± 0.005, Im(a) = 0.00 ± 0.08, 
c = 0.018 ± 0.008 and d = —0.059 ± 0.012, where the errors are statistical only. Although 



the fitted results are consistent with world average values 18| , possible correlations between 



the X and Y are not considered. Also the fitted value of the parameter d is not consistent 
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FIG. 3: (a) The experimental Dalitz diagram for the decay r]' — )• r]TT~^7r~ in terms of the variables 
X and Y with the r]Tr^TT^ mass in the ij' mass region. The corresponding projections on variables 
X and Y are shown in (b) and (c), respectively, where the dashed histograms are from MC signal 
sample with rj' — t- rjn^Ti^ events produced with phase space and the solid histograms are the fitted 
results described in the text. 

with zero so the matrix element can not be well described by a linear function of Y only. So 
we do the fits to the Dalitz plot described below. 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 




-1 



-0.5 




X 



0.5 



1 



0.35 

0.3 I 
0.25 r 

0.2 r 
■0.15r 

0.1 r 
0.05 r 

0-1 




-0.5 







0.5 



1 



FIG. 4: Dependence of the square of the r/' decay matrix element on the Dalitz variables X and 
Y . The solid lines are the results of the fits of the data described in the text. 



X [N, a, 6, c,d) = 2 



(4) 



In the fitting procedure, the Dalitz plot is subdivided into 26 X-bins and 22 F-bins, i.e. 
572 cells in total. Dalitz plot parameters are obtained by minimization of the function: 

■I ' 

Here the index i enumerates cells in Dalitz plot (empty cells outside the Dalitz plot bound- 
aries are excluded), N is normalization factor, a,b,c and d are the Dalitz plot parameters. 
The Mi and Di are the numbers of (weighted) entries in the i-th bin of the two-dimensional 
histograms in the Dalitz variables for MC and for the background-subtracted data, respec- 
tively. The statistical error a includes background subtraction and MC statistical errors. 
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The MC histogram is obtained as follows: 

Nev 

Mi = ^(1 + aYj + hY^ + cXj + rfXj), (5) 
i=i 

for the general decomposition parametrization, where the index j is over the generated 
events and Xj and Yj are the generated true values of Dalitz variables. Similarly for the 
linear parametrization, 

Nev 

M, = "^{\l + aYj\^ + cXj+dX]). (6) 
i=i 

The fit procedure has been verified with MC by checking the input and output values of the 
Dalitz plot parameters. 

First we fit using the general decomposition parametrization of the matrix element and 
obtain the following values for the parameters of the matrix element and for the correlation 
matrix {x^/NDF = 504/476, where NDF is the number of degrees of freedom.): 



a = -0.047 ±0.011 



/l.OOO -0.442 -0.010 -0.239\ 



b= -0.069 ±0.019 
c = +0.019 ±0.011 
d= -0.073 ±0.012 



(7) 



1.000 0.025 0.282 
1.000 0.030 
\ 1.000 / 

The errors are statistical only. This result is illustrated in Fig. [5l where we show the com- 
parison of data (dots with error bars) and MC weighted with fitted coefficients (histogram) 
as a function of Y, in different X- intervals for rj' — t- rj-K^-K^ . Parameter c is consistent with 
zero within I.80". The fitted results are almost the same with the value of parameter c fixed 
at zero. The statistical significance of c is estimated to be 2. la, from the difference of the 
value taking the difference in the number of degrees of freedom (ANDF = 1) in the fits 
into account. 

The extra term eXY or fX^ + gY^ has also been added into the general parametrization. 
The fitted value of parameter e is 0.000 ±0.018, which is consistent with the conclusion from 
the VES measurement [9|. The fitted results of parameters f and g are 0.037 ± 0.035 and 
—0.014 ± 0.018, respectively, and the corresponding statistical significances are very small 
(~ la). All the other parameter values are almost the same. 

We also perform a fit using the linear parametrization of the matrix element and obtain 



13 



250 7 
200- 
150r 
lOOr 
50 r 

O 250- 



(0 



200 r 



§ 150r 
50 r 



1.1<IXI<1.2 



0.7<IXI<0.8 



250 r 

200^ 

150r 

100- 

50 r 





,0.3<IXI<0.4 
II I I I 



1.0<IXI<1.1 



0.6<IXI<0.7 



0.2<IXI<0.3 
I I 



,0.9<IXI<1.0 



0.5<IXI<0.6 




,, ,0.1<IXI<0.2 

I I 



0.8<IXI<0.9 



0.4<IXI<0.5 

iMl'l 



IXI<0.1 



-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 



FIG. 5: Experimental distributions of the variable Y in various intervals of X with the fitting 
function (histogram) for the general decomposition parametrization. 



(xVNDF = 521/476): 
Re{a) 



-0.033 ±0.005 



/l.OOO -0.001 0.001 -0.138\ 



Im(a) = 0.000 ± 0.049 
c = +0.018 ±0.009 
d= -0.059 ±0.012 



(8) 



1.000 0.000 0.000 
1.000 0.024 
\ 1.000 / 

The errors are statistical only. The parameter c is consistent with zero within 2.0cr, and the 
statistical significance is estimated to be 2.2a. 



VII. SYSTEMATIC UNCERTAINTIES 

The sources of the systematic errors for the branching fraction measurement are summa- 
rized in Table [B The uncertainty is negligible for pion identification since the identification 
of only one of the pions is required. The uncertainty in the tracking efficiency is 1% per 
track and is additive. The uncertainty associated with the kinematic fit is determined to be 
0.2% using the control sample J/ip ^ n^n^n^. The uncertainty due to photon detection 
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is 1% per photon. This is determined from studies of photon detection efficiencies in well 
understood decays such as J/tp ^ p^vr^ and study of photon conversion via e+e^ — t- 77. Ac- 
cording to the MC simulation, the trigger efficiency for signal events is almost 100%, and the 
uncertainty is neglected. The background uncertainties are evaluated by changing the back- 
ground fitting function from a first order polynomial to second order and the fitting range, 
resulting in change of the branching fraction by 0.3%. The uncertainties of B{ri' rjTc^TT') 
and B{ri — )■ 77) are 1.6% and 0.5%, respectively |18!]. The fitted results to the Dalitz plot 
matrix element show correlations between the Dalitz plot parameters. This should be prop- 
erly taken into account when integrating the amplitude over phase space to obtain the decay 
width. The maximum difference in the efficiency is 2.5% by using general decomposition 
parametrization results or linear parametrization results. The difference (2.5%) is conserva- 
tively taken into the systematic errors. Finally the uncertainty on the number of J/ip events 
is 2%. Assuming that all of these systematic error sources are independent, we add them in 
quadrature to obtain the total systematic error shown in Table HI 

TABLE I: Relative systematic errors (%) for the branching fraction measurement. 



Somxe 


B{J/i; ^ W) 


Part ID 




Tracking 


2.0 


Kinematic fit 


0.2 


Photon efficiency 


3.0 


MC statistics 


0.3 


Trigger efficiency 




Background shape 


0.3 


Intermediate branching fractions 


1.7 


DaUtz plot matrix element 


2.5 


Number of J ftp events 


1.3 


Sum in quadrature 


4.9 



The systematic errors in the measurement of the Dalitz plot matrix element are summa- 
rized in Table mi The uncertainty from the backgrounds is negligible since the contamination 
is very small (~0.57%). The tracking efficiency correction functions for tt~^ and tt~ are ob- 
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tained by using the control sample J /ip ^ 7[~^7T~pp, where the transverse momentum region 
of pion has covered the region of signal pion transverse momentum. The differences on 
the fitted values of parameters a, b, c and d are 3.3%, 3.3%, 4.4% and 1.2% in the general 
parametrization, and 3.3%, 4.8%, 2.4% for the parameters Iie{a), c and d in the linear 
parametrization by applying the tracking efficiency correction functions for 7r+ and tt~, re- 
spectively. The differences on the fitted results of parameters due to changing the tjit'^it' 
mass requirement are included in the systematic errors. The fitted results are also compared 
using a 4C instead of the 6C kinematic fit, and the corresponding differences are taken as 
the systematic errors due to the kinematic fit uncertainty. Binning size was changed up 
to a factor of two: 0.1 < AX, AF < 0.2. The biggest differences on the fitted parameter 
values are taken as the systematic errors due to the binning size uncertainty. To determine 
the systematic errors associated with the event selection, especially for the selection of r) 
candidates, another set of event selection criteria are applied: (1) The photon with the max- 
imum energy is regarded as the radiative photon {'jrad)- (2) We do a 4C kinematic fit to the 
J/ip ^ 7rad777r^7r~ hypothesis looping over all the other photon candidates and select the 
combination with the minimum X^(7rad777r^7r^). (3) The 77 mass is required to be within 
the rj mass region. After applying the above event selection criteria, the difference on the 
total number of signal events is only about 0.53%. The fits to the Dalitz plot parameters 
are done with these events, and the differences are included into the systematic errors due 
to the event selection method uncertainty. The PHOTOS package |23| was used to include 
final state radiation (FSR). By changing the ratio of FSR events, the differences are taken 
into the systematic errors due to the FSR simulation uncertainty. Assuming that all the 
sources are independent and adding them in quadrature, one gets the total systematic errors 
of parameters a, b, c and d in the generalized representation are 4.9%, 12%, 12% and 3.1%, 
and the total systematic errors of parameters Re(a), c and d are 7.2%, 12% and 6.5% in the 
linear representation, respectively. 
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TABLE II: Relative errors of the parameters of the matrix element for the generalized and linear 
represent at ions . 



Source 


generalized representation 


linear 


representation 




Qj 


u 


c 


rl 


Re(a^ 


c 


rl 


Tracking efficiency 


3.3 


3.3 


4.4 


1.2 


3.3 


4.8 


2.4 


m„^+^- mass cut 

I/TV ' TV 


0.9 


4.8 


3.3 


1.7 


1.4 


2.1 


3.2 


Kinematic fit 


2.8 


4.9 


2.1 


0.7 


5.2 


7.0 


4.7 


Binning size 


0.9 


8.0 


9.2 


1.5 


2.6 


5.7 


1.7 


Different selection method 


1.6 


2.9 


0.7 


1.4 


2.0 


4.8 


1.1 


FSR simulation 


1.0 


0.4 


1.9 


0.2 


0.9 


0.7 


0.4 


Sum in quadrature 


4.9 


12 


12 


3.1 


7.2 


12 


6.5 



VIII. SUMMARY 



Using the large J/ijj sample ((225.2 ± 2.8) x 10^ J/ip events) collected with BESIII, the 
branching fraction of J/ip 777' is measured to be 

B{J/ip 77/') = (4.84 ± 0.03 (stat) ± 0.24 (sys)) x 10-^ 



which is consistent with the recent BESII value ((5.55 ± 0.44) x 10 ^) 2J] within 1.5cr, and 
the CLEO value ((5.24 ± 0.17) x 10~^) [25I] within 1.4cr, which are used in obtaining the 
world average value (5.28 ± 0.15) x 10^^) by the PDG |l8|. 

The parameters of the matrix element for the decay process 77' — t- t^tt+tt^ have been 
determined for the generalized and linear representations. They are: 



a = -0.047 ±0.011 ±0.003 
b= -0.069 ±0.019 ±0.009 
c = ±0.019 ±0.011 ±0.003 
d = -0.073 ±0.012 ±0.003 
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for the generalized parametrization, and 

Re{a) = -0.033 ± 0.005 ± 0.003 
Im(a) = 0.000 ± 0.049 ± 0.001 
c = +0.018 ±0.009 ±0.003 
d = -0.059 ±0.012 ±0.004 

for the hnear parametrization, where the first errors are statistical and the second systematic. 

Table IIIII shows the experimental and theoretical values of the parameters of the matrix 
element squared for t]' rj-K^ix^ in the general parametrization (second, third and fourth 
columns) and in the linear parametrization (sixth, seventh and eighth columns). The the- 
oretical values in Ref. [2fl] are the latest calculations within the framework of f/(3) chiral 
effective field theory in combination with a relativistic coupled-channels approach. We see: 
(1) The errors of our fitted parameter values are smaller compared to previous published 
results. (2) In the general decomposition parametrization of the matrix element, the central 
values of parameters a and h are consistent with the results from GAMS-47r Collabora- 
tion [10], where the neutral decay rj' rjii^TT^ events were analyzed; however the central 
values of parameters c and d are consistent with the results from VES Collaboration [9]. (3) 
The negative value of the coefficient b indicates that the two kinds of parametrization are 



not equivalent. This conclusion is consistent with that from GAMS-47r Collaboration 10| ; 



however it is different from the conclusion from the VES Collaboration [9|, where the fit 
with linear parametrization yields an unsatisfactory /NDF = 170.5/114 ratio. (4) The 
quadratic term in X is unambiguously different from zero. Similarly for the quadratic term 
in Y. The measured value of the F-variable quadratic term (b) is not consistent with the 
expected value of around zero in the Effective Chiral Lagrangian model in which the lowest 
lying scalar meson candidates a and k together with the /o(980) and ao(980) are combined 
into a possible nonet 28j; however it can be accommodated in a U{3) Chiral Unitarized 
model by including final state interactions j3|. The dynamical nature of this term needs fur- 
ther clarification. (5) The value of parameter c, which tests C parity violation in the strong 
interaction, is consistent with zero within 2a in both parametrizations. In the future, with 
much more BESlll data, the hadronic decays of rj' can be measured with higher precision, 
especially the Dalitz decay parameters, allowing more stringent testing of the predictions of 
ChPT [291. 
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TABLE III: Experimental and theoretical values of the parameters of the matrix element squared 
for rl — )• iqii^-K^ in the general parametrization (second, third and fourth columns) and in the 
linear parametrization (sixth, seventh and eighth columns). 



Par. 


VES [9] 


Theory [26j 


This work 


Par. 


CLEO [nj 


VES [27] 


This work 


a 


-0.127 ±0.018 


-0.116 + 0.011 


-0.047± 0.012 


Re(a) 


-0.021 + 0.025 


-0.072 + 0.014 


-0.033+0.006 


b 


-0.106 ±0.032 


-0.042 ±0.034 


-0.069 ±0.021 


Im(a) 


0.000 (fixed) 


0.000+0.100 


0.000 + 0.050 


c 


+0.015 ±0.018 




+0.019+0.012 


c 


0.000 (fixed) 


+0.020+0.019 


+0.018+0.010 


d 


-0.082 ±0.019 +0.010 ±0.019 


-0.073 + 0.013 


d 


0.000 (fixed) 


-0.066 + 0.034 


-0.059+ 0.013 
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